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ABSTRACT 

The random-anisotropy Blume-Emery-Griffiths model, which has been pro- 
posed to describe the critical behavior of '^He-'^He mixtures in a porous medium, 
is studied in the pair approximation of the cluster variation method extended to 
disordered systems. Several new features, with respect to mean field theory, are 
found, including a rich ground state, a nonzero percolation threshold, a reentrant 
coexistence curve and a miscibility gap on the high ^He concentration side down 
to zero temperature. Furthermore, nearest neighbor correlations are introduced in 
the random distribution of the anisotropy, which are shown to be responsible for 
the raising of the critical temperature with respect to the pure and uncorrelated 
random cases and contribute to the detachment of the coexistence curve from the 
A line. 
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1. Introduction 

The study of the effects of randomness on phase transitions has a long 
story and only recently there has been considerable experimental and theo- 
retical efforts to deeply understand them. Here we will be concerned with a 
simple model for the phase separation of ^He-^He mixtures in aerogel[l]. 

Aerogel is an extremely porous medium, the silica glass, made via the 
sol-gel process[2]. Its porosity may be as high as 99.8%[3] corresponding to 
a density of 4mg/cm^. There is experimental evidence[2,4] that the aerogel 
microstructure is rather ramified, composed of silica strands with a thickness 
of order 10 A and with an area of order 700m^/gr[5]. 

In order to describe in a simple way the effects of the porous medium on 
the phase separation diagram of ^He-^He we model the aerogel as a random 
external field which selects which of the two types of Helium to have nearby. 
Thus it is quite interesting, among various other things, to understand the 
role of the two types of impurities, i.e. the annealed ^He and the quenched 
randomness of the external field, on the superfiuid transition of ^He. 

Using mean field theory, transfer matrix and real space renormalization 
group calculations it was possible to predict a variety of physically acceptable 
scenarios for the phase separation diagram as the density of the quenched im- 
purities is varied[l]. The main features are the disappearance of the tricritical 
point present only in the pure ^He-'^He system[6] and superfiuidity in two co- 
existing phases, one rich in ^He and the other in ^He. Recent experiments [7] 
have confirmed the above qualitative picture. 

The goal of this paper is to give a detailed study of the model proposed in 
[1] using the pair approximation of the cluster variation method (CVM)[8], 
generalized to random systems [9, 10]. This approximation improves the stan- 
dard mean field by taking into account effects of correlation among spins 
at different sites. Most of the qualitative features of the mean field phase 
diagram are preserved. 

Our approach allows us to incorporate effects of nearest neighbors (n.n.) 
correlations in the quenched impurities, which are essential to a full under- 
standing of the role of aerogel in the •^He-'^He phase diagram. Unlike the 
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coexistence curve, the A line separating the normal fluid from the superfluid, 
is rather sensitive to correlations of the quenched impurities at high ^He con- 
centration. Interesting new features are found in the T = phase diagram 
and the role of the percolation threshold is elucidated. 

In the next section the random-anisotropy Blume-Emery-Griffiths model 
(RABEG) will be introduced. The CVM generalized to random systems will 
be discussed in detail in connection with the RABEG model. Sec. 3 is devoted 
to the discussion of the phase diagram obtained in the present approximation, 
also in connection with experimental results, and in Sec. 4 some conclusions 
are drawn. 

2. Definition of the model and the CVM pair approximation 

The bulk phase diagram of ^He-^He is well described by the Blume- 
Emery-Griffiths (BEG) model[6] which is a spin-1 system with nearest neigh- 
bor interactions and a uniform anisotropy external field distinguishing only 
between the ±1 and values of the spin. Its generalization to the ^He-^He 
mixture in a porous medium is described by the hamiltonian[l] 

h = -jY^ s,s, -kY^ sfs] + J2 ^^sf, (1) 

where the spin variables Si = 0, ±1 are defined on sites of a lattice A. The 
first and the second sums are over n.n. pairs. When dealing with ^He-^He 
mixtures[6] Si = ±1 and = represents ^He and ^He respectively, and 
the superfluid transition is associated to the spontaneous symmetry breaking 
between the ±1 states. The exchange interaction with strength J > is 
responsible for the superfluid ordering, while the random anisotropy Aj is 
related to the difference of chemical potentials fj,s — fj,4 at site i, and K = 
K33 + K44 — 2^34, where Ka/s is the interaction energy between "He-''^He 
atoms. Koii3 does not depend significantly on a and /? and then it is generally 
assumed K = 0. 

Since we are going to study a model for ^He-^He mixtures in aerogel, we 
will let each Aj take the value Aq or Ai with the same meaning as in [1], i.e. 
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Ao < is the value of the anisotropy at the pore-grain interface of aerogel 
(where ^He prefers to stay), while Ai can be thought of as a bulk field which 
controls the ^He concentration x. 

Two order parameters can be defined: m = {{Si))A and q = {{Sf))A, 
which are the quenched averages of the thermal equilibrium values (Si) and 
(Sf). m is the superfluid order parameter, while q = 1 — x is the ^He 
concentration. 

As we mentioned in the introduction we take into account n.n. correla- 
tions in the random field distribution assuming the following joint probability 
density (with i and j nearest neighbors) for the Aj's: 

P(Ai,A,)= Yl Pn,nJ{Ai-AnM^j-^n,), (2) 

"1,^2=0,1 

Pnma = PmPn2 + (2(^nin2 " 1)^/, (3) 

where po = p is the fraction of sites at the aerogel interface (the grain space 
is neglected), pi = 1 — p and e is a measure of the interface correlation. 
Integrating Aj one obtains 

P(A,) = pdiA, - Ao) + (1 - p)d{A, - Ai), (4) 

which was used as a starting point in [1] . 

We are now ready to build the pair CVM free energy following the pro- 
cedure outlined by Morita[9] for a general random system, which we briefly 
review. 

Consider a random system on a finite lattice, for which the configuration 
of random fields and interactions is specified by a unique random variable h 
with distribution P{h), and define, according to the usual rules of statistical 
mechanics, a density matrix pc{a\h) for each configuration h {a stands for 
the set of statistical degrees of freedom of the system, e.g. spin variables). 
In terms of pc, the free energy associated to the configuration h will be given 

by 

F{h) = J2pc{<j\h) [H{h, a) + kBTlnp,{a\h)] , (5) 
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where H{h, cr) is the hamiltonian, ks Boltzmann's constant, and T absolute 
temperature. The quenched free energy F will then be given by 

F = J2P{h)F{h)- (6) 

h 

Upon introducing the generalized density matrix 

p{h, a) = P{h)pc{a\h) 

one can easily show that 

F = Y, P{h, (^) [H{h, a) + ksT In p{h, a)] + T^c, (8) 

where Sc depends only on P{h) and is given by 

h 

As in the pure case, we have a variational principle: p(/i, cr) can be determined 
as the matrix which minimizes F, as given by (8), with the constraint 

Y,P{KcT) = P{h). (10) 
cr 

In such a scheme it can be shown that the quenched average of the expectation 
value {A{h^ a)) of an operator A{h^ a) is given by 

{{A{h,a)))h = J2Pih){A{h,a)) = Pi^, a)A{h, a). (11) 

h h,a 

The cluster variation method can then be obtained by taking the ther- 
modynamic limit and truncating the cumulant expansion for the entropy 
S = — ksT'^f^ ^ p{h, a) In p{h, a) to a set of "maximal preserved clusters" 
Fj, i = 1,2, ...r (and all their translates). The variational principle will 
then be applied to the reduced density matrices pr^ {h, cr) associated to the 
maximal preserved clusters. 

In the following, our maximal preserved cluster will be the n.n. pair, 
and thus we will introduce a density matrix pp with elements pp{h,a) = 

5 
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p^p"'''\s^,S2) (ni,n2 = 0,1, SuS2 = ±1,0), where h = (Ani,AnJ is the 
random field configuration on the cluster and a = {81,82) denotes the set 
of spin variables, pp (which is diagonal, since we are dealing only with the 
z-component of the spins) is subject to the constraint "^a- Pp{h, a) — P{h), 
that is 

J]p(--"^)(5i,52)=p„,„,. (12) 
S1S2 

Writing direct sum of density matrices, one for each possible config- 

uration of the random fields on the cluster: 

Pp= Pn.n.P;"-"^) (13) 
ni ,n2=0,l 

the constraint (12) becomes 

IVppK'-^) = 1, ni,n2 = 0,l. (14) 

Furthermore, pp has the obvious symmetry property 

pj,-i'-^)(5i,52) = pJ."^'"^n'52,-Si) > 0. (15) 

As usual, a (reduced) site density matrix can be obtained from pp by a partial 
trace: 



or, more explicitly. 



Ps= Tr Pp, (16) 

p\s 



p(-i)(^i)= ^p(-i.-^)(5i,52). (17) 

n2S2 



The pair CVM free energy is then given by 



{l-u)Tr{pAnps)+'^Tr{pplnpp) , (18) 



with U = {{H))/^ and v coordination number of the lattice, and has to be 
minimized with respect to Pp, with the constraints (14) and (15). Taking 
derivatives of F with respect to pp leads to three systems of equations, which 
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are linear in the elements of lnpp^°°\ Inpp^^"'^'' and Inpp*^^^^ respectively, and 
where the elements of In pg can be treated as parameters. Exponentiating 
the solution and making use of (14) yields 



n2 



81,82 



(19) 



with (3 — lyJ/ksT and 



= eu' 



2{l-qn) 



Vr 



(n) 



where a = {1^ — 1) /i/, dn = A„/ (lyJ). The parameters and qn, which can 
be thought of as "local order parameters" at the aerogel interface (n = 0) 
and in the interior (n = 1), are given by 



9n = p1"n+)+Pl"n-) 



(20) 



It can be easily shown that the parameters m„ and q„ are solutions of the 
following equations 



"2,^2 

Pnqn= J2 (p?""^^( + ,'52)+p("''^^)(-,'52)) 



(21) 



n2,S2 



which, together with (19), are the basic results of this section, and the starting 
point for the analysis of the phase diagram of our model at finite temperature. 

Finally, the order parameters of the model are related to our " local order 
parameters" by 

m = pmo + (1 — p)mi 
q = pqo + il-p)qi 



(22) 
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3. The phase diagram 

At zero temperature, several different ptiases can be found, wliicli will be 
identified by the set (moo = 9oO)"^oi = 9oi)"7,io = qiOjiTin = qn), where 

S 
S 

These are related to the expectation values of Si and Sf respectively when 
the anisotropy field is at site i and A^^ in one of its neighbours. 

The phase diagram is easily determined by comparing the energies U of 
the different phases, and is reported in Fig. 1. Two comments are in order: i) 
as opposite to the mean field result [1], the boundaries between the different 
ground states do no longer depend on the probability distribution of the 
random variables and ii) increasing the level of the approximation we have 
obtained a more complicated phase diagram, with a larger number of possible 
ground states. This is reminescent of what happens in the exact solution at 
T = of the random field Ising model on a Bethe lattice[ll]. Notice that our 
model (1) reduces to a random field Ising model for J = and K ^ since 
all spin variables appear squared and thus assuming only the two values 0, 1. 

At finite temperature, and for Aq < 0, it has been shown in the mean 
field analysis [1] that the phase diagram has a quite rich structure, with a 
second order transition separating the ordered phases from the disordered 
one, a first order transition, terminating in a critical point, between the 
ordered phases above, and several multicritical points. Our approximation 
yields an even richer structure: indeed we have three different ground states 
(see dotted line in Fig. 1), (1111) for di < 1/2 corresponding to ^He present 
everywhere, (1110) for 1/2 < di < 1, which should correspond to ^He present 
in the aerogel interface and diffusing near it, and (1100) for di > 1, i.e. ^He 
only at the aerogel interface. Again, ordered phases are separated by the 
disordered (high temperature) one by a second order transition line, for which 
an equation can be obtained by expanding (23) around mo = mi = 0. Up to 
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(23) 



the first order, the equations for yield 

mo = aooniQ + aoimi 
777-1 = ciiomi + aiimi 

where 



5^A + 4^^ sinh^ 

p^qx Z^x i^. 



(24) 



(25) 



The equation for the critical temperature turns out to be 

(aoo - l)(aii - 1) = aoiaio, (26) 

meaning that a non trivial solution of eq. (24) exists, where go and qi must be 
determined from the corresponding equations with uiq = mi = 0. A solution 
for the critical temperature is obtained for any di only if p is greater than 
the percolation threshold Pc- 

The percolation threshold can be defined as follows: in the limit Ai — > 
+00, Si will be for all sites with anisotropy Ai and our model will be equiv- 
alent to a random-diluted BEG model with anisotropy Aq and concentration 
p. Such a model will undergo a second order transition at Tc > for any 
large enough p, and the percolation threshold is just that value of p for which 
Tc becomes zero. It can be obtained by taking the limits Ai +oo and 
P — > +00 in the equation for the critical temperature, and the result is 

representing a remarkable improvement with respect to the mean field theory, 
which gives pc = 0. Notice that for uncorrelated disorder, i.e. e = 0, the exact 
result for the Bethe lattice with coordination u is recovered [12]. 

For p > Pc the phase diagram has the structure shown in Figs. 2(a) and 
2(b) (heavy lines): solid lines are used for the second order transition, while 
dashed lines stand for the first order transition separating (1111) and (1110) 
phases (if in the {di,T = plane) or for the corresponding coexistence 

region (if in the (x, r) plane, where x is the ^He concentration). C is a critical 
point. No first order line separating (1110) and (1100) phases seems to start 
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from the point F{di = 1,t = 0), where a first order transition occurs at 
zero temperature. Again this is reminescent of the behaviour of the random 
field Ising model on the Bethe lattice[ll]. The phase diagram at T > is 
qualitatively the same as that given by mean field theory, which is reported 
for comparison (light lines). In all considered cases (Aq < 0) x < 1 — p. 

For < p < pc Si. new feature arises, that cannot be described by mean 
field theory: the second order line has a limit point F at zero temperature 
and (ii = 1, independent of p, e and do, and the critical temperature is zero 
for any di > 1, meaning that, even if the ground state is ordered, aerogel 
concentration is too low for the pore-grain interface to sustain an ordered 
phase at finite (i.e. nonzero) temperature. The corresponding phase diagrams 
are given in Figs. 3(a) and 3(b) and in Figs. 4(a) and 4(b), for typical values 
of the parameters. In the latter case the second order line is reentrant and 
intersects the first order line in the critical end points A and B. It is easily 
realized that there exists a value p* < Pc oi p for which the points A and 
B merge into one, giving rise to a new multicritical point where the second 
and the first order lines are tangent with respect to each other. It is also 
interesting to remark that the first order transition has always a reentrant 
behavior, a feature that does not exist in the mean field treatment [1]. 

In our scheme the behavior of order parameters and n.n. correlation 
functions can be obtained in a natural way: in Fig. 5 we report the behavior of 
the local order parameters rrin and vs. temperature for the case 1/2 < di < 
1 (dotted line in Fig. 3(a)), where mo = qo = 1 and mi = qi = Pio/(l — p) 
in the ground state. It is interesting that mo and mi exhibit a quite unusual 
oscillation. Furthermore, in Fig. 6, we report the global order parameter m 
vs. temperature at fixed concentration, corresponding to the dotted line in 
Fig. 2(b). 

Finally, let us discuss the effect of a correlation e > in the random distri- 
bution on the phase diagram, and especially on the second order transition. 
In the mean field analysis[l] it was shown that the second order transition 
occurs at T = 1 — x = independent of p, while our analysis shows that 
the second order critical temperature at fixed concentration is indeed weakly 
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dependent on p (Fig. 7). Introducing a correlation e of order unity in the 
random distribution of the Aj's causes the critical temperature to raise sig- 
nificantly, especially in the \ow-q (high-x) region (Fig. 8), and we believe 
that this effect should account for the experimentally observed [7] increase in 
the critical temperature with respect to the pure case. The behavior of Tc as 
a function of e is shown in Fig. 9. 

4. Conclusions 

A model for the critical behavior of ^He-^He mixtures in a porous medium 
(aerogel), namely the random-anisotropy Blume-Emery-Griffiths model[l], 
has been investigated in the pair approximation of the cluster variation 
method. This improves the mean field analysis given in [1] in several di- 
rections. A rather rich phase diagram is found already at zero temperature, 
with various ground states corresponding to (Aq < 0): ^He present every- 
where, (1111) in Fig. 1, at the aerogel interface and in a neighbour of it 
(1110), and only at the aerogel interface (1100). The symmetric situations, 
with ^He preferentially near aerogel are of course possible within the context 
of a model (Fig. 1 for do > 1/2). 

A nonzero percolation threshold can be determined, above which the 
phase diagram is mean field like, and below which it changes qualitatively, 
with the appearance of a zero temperature limit point for the second order 
transition line. A completely open and extremely interesting question con- 
cerns the universality class of the whole line of continuous transitions: is 
there a zero temperature or di = +oo fixed point attracting a finite part of 
the critical line and an unstable fixed point at finite temperature separat- 
ing two different critical behaviors (e.g. di < and di 1)? A real space 
renormalization group analysis might be suitable to address this problem. 

In all considered cases, the first order transition line (in the temperature- 
anisotropy phase diagram) and the coexistence region (in the temperature- 
concentration phase diagram) exhibit reentrant behavior. This reentrance is 
also partially responsible for the existence of a miscibility gap on the high 
^He concentration side going down to zero temperature. 
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The experimental phase diagrams are in remarkable agreement with the 
ones presented in figs. 2(b) and 3(b). Quantitative comparisons would of 
course require more adequate treatments of the superfluid transition. 

Particular attention has been devoted to the analysis of the effects of a 
nearest neighbor correlation in the random distribution of the anisotropy, 
which should not be neglected due to the high correlation characterizing the 
aerogel interface. The main result is the increase in the critical temperature 
with respect to the pure case, contributing to the detachment of the A line 
from the coexistence curve and which has been experimentally observed [7] 
and cannot be explained in terms of the uncorrelated randomness alone. 
Correlations in the quenched impurities make the question of the universality 
class of the critical line more involved even at low disorder concentration. 
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Figure captions 



Fig. 1 : Phase diagram at T = 0. 

Fig. 2(a) : Phase diagram in the (cii, r) plane for = 6, e = 0, p = 0.21 > 

do = —0.5. Heavy hnes are our results and light lines are results 
from mean field approximation. 

Phase diagram in the (x, r) plane for the case of Fig. 2(a). Sym- 
bols as in fig. 2(a) (for the dotted line see Fig. 6). 
The same as Fig. 2(a) with p — 0.19 < Pc (for the dotted line see 
Fig. 5). 

Phase diagram in the (x, r) plane for the case of Fig. 3(a). 
The same as Fig. 2(a) with p — 0.1 < Pc- 
Phase diagram in the {x,t) plane for the case of Fig. 4(a). 
rrin and Qn vs. T along the dotted line in Fig. 3(a) {di = 0.65). 
m vs. T along the dotted line in Fig. 2(b) {x = 0.69). 
The critical temperature in the (x, r) plane for = 6, e = 0, 
do = —0.5 and p = (solid line), 0.1 (dashed line) 0.2 (dotted 
line) . 

The critical temperature in the (x, r) plane for v = 6, do = —0.5, 
p = 0.1 and e = (solid line), 0.5 (dotted line), 1.0 (dot-dashed 
line) and 2.0 (dashed line). 
Fig. 9 : The critical temperature vs. e for the case of Fig. 7 and with 
X = 0.6. 
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